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The flow instability and further transition to turbulence in a toroidal pipe (torus) with 
curvature (tube-to-coiling diameter) 0.049 is investigated experimentally. The flow inside 
the toroidal pipe is driven by a steel sphere fitted to the inner pipe diameter. The sphere 
is moved with constant azimuthal velocity from outside the torus by a moving mag¬ 
net. The experiment is designed to investigate curved pipe flow by optical measurement 
techniques. Using stereoscopic particle image velocimetry, laser Doppler velocimetry and 
pressure drop measurements the flow is measured for Reynolds numbers ranging from 
1000 to 15 000. Time- and space-resolved velocity fields are obtained and analyzed. The 
steady axisymmetric basic flow is strongly influenced by centrifugal effects. On an in¬ 
crease of the Reynolds number we find a sequence of bifurcations. For Re = 4075 ± 2% a 
supercritical bifurcation to an oscillatory flow is found in which waves travel in streamwise 
direction with a phase velocity slightly faster than the mean flow. The oscillatory flow 
is superseded by a presumably quasi-periodic flow at a further increase of the Reynolds 
number before turbulence sets in. The results are found to be compatible, in general, 
with earlier experimental and numerical investigations on transition to turbulence in 
helical and curved pipes. However, important aspects of the bifurcation scenario differ 
considerably. 


1. Introduction 

Transition to turbulence in straight circular pipes is one of the oldest and most fun¬ 
damental problems of flu id mechanics, e ntailing various fundamental questions about 
the nature of turbulence ( Eckhard^l2008ll . The understanding of the physics of transi¬ 
tion to turbulence in straight pipes has experienced a significant progress in recent years 
due to the application of mo dern optical measurement techniques and computati onally- 
based theoretical modeling ( Hof et all 2004 : Eckhardt et all EoOTI : Mnllinl [201 ih . It is 


characteristic for straight pipe flow that transition occurs despite of the linear stability 
of the laminar Hagen-Poiseuille flow if the flow perturbations exceed a certain thresh¬ 
old. Moreover, and different from plane Poiseuille flow, Rayleigh-Benard convection, or 
Taylor-Couette flow, there exists no instability of Hagen-Poiseuille flow for any Reynolds 
number Re = Ud/v where U is the mean velocity, d the diameter of the tube and v the 
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kinematic viscosity of the fluid ( Drazin fc R.eidI 198l[) . Cons equently, turbule nce can only 
be triggered by flnite amplitude perturbations. Following iR.evnol^ ( 188,‘lh the critical 
point is then defined as the Reynolds number above which turbulence is first sustained 
indefinite ly, whereas below any turbulence initially in the flow will eventually decay. As 
shown bv lAvila et al. ( 201lll this critical Reynolds number is RCc ~ 2040 where the decay 
of turbulence is balanced by a spreading process. 

In curved pipes the process of transition to turbulence qualitatively differs from that in 
straight pipes. Curved-pipe flow has received less attention than straight-pipe flow, even 
though the former is of no lesser importance: practically all pipes in nature, e.g. blood 
vessels, or in engineering are curved. In curved pipes the basic flow is strongly affected 
by an imbalance between the cross-stream pressure gradient and the centrifugal force 
which leads to a secondary cross-stream motion, typically in form of a pair of steady 
streamwise Dean vortices symmetric with respect to the tangent plane. With increasing 
Dean number the maximum of the streamwise velocity is shifted radially towards the 
outer wall of the pipe. Owing to the increased gradient of the streamwise flow the drag 
in curved pipes is considerably hig her than in straight p ipes. For reviews of curved-pip e 
flow and further reading we r e fer to Berger et al\ ( 1983h : lHuttl fc Friedrich ( 20001 [200lh : 
Naphon fc Wongwised ( 2006l l: IVashisth et al (200^! 


Dean ( I927I . 1928 1 has solved the simplified Navier-Stokes equations for a coiled pipe 


of small curvature showing that the flow is governed by two parameters. These param¬ 
eters are the curvature S = d/ D of the pip e diameter d to the coiling diam eter D and 
the Dean number De = R,e\/6. iTavlor ( 1929f l. lwhitel (192^ and Adleil ( 1933 l found that 
the flow in curved pipes remains laminar up to Reynolds numbers higher at least by a 
factor of two than in straight pipes. They also noticed that transition to turbulence is 
not as abrupt as in straight pipes but occurs graduall y, without any discontinuity of the 
characteristic observables such as the pressure drop. Sreenivasan fc Strvkowski ( 1983ll 
investigated curved-pipe flow for moderate Reynolds numbers. They found that the tur¬ 
bulent flow in a straight pipe becomes laminar after entering a coiled-pipe of the same 
diameter. They recognized, moreover, that the curved-pipe flow becomes time periodic 
before turbulence commences. The oscillation amplitude was found to be strongest near 
t he inside wall of the pip e. _ 

Webster fc Humphrey ( I993ll investigated the unsteady three-dimensional flow through 


a helical pipe at transitional conditions using laser-Doppler velocimetry. For 5 = 0.055 
they measured the onset of periodic low-frequency perturbation waves at Rcc = 5060. 
They specified the nondimensional frequencies found with fd/U = 0.25 and 0.5 for all 
Re ynolds numbers investi g ated. Complementary numerical simulations for Re = 5480 
by Webster fc Humphre y! (Il99it) were in agreement with the experimental results of 
I Webster fc Hump hrey (1993 ). Dye str e aks, u sed for visualization, were found to diffuse 
at Re = 6330. Webster fc Humphrey ( 1997t) interpreted this behavior as the onset of 
turbulent fluctuations. 

Flows in weakly curved ducts with constant curvature and in helical pipes with small 
torsion have a common asymptot ic limit: toroidal- p i pe flo w. However, the inlet and outlet 
conditions are generally different. Piazza fc Ciofalol ( 20Ilh numerically simulated the flow 
in a circular toroid with a circular cross-section which is driven by a constant streamwise 
body force, using periodic boundary conditions. For 5 = 0.3 a supercritical Hopf bifurca¬ 
tion was found in the interval 4556 < Red < 4605, giving rise to an azimuthally traveling 
wave which took the form of a varicose streamwise modulation of the two Dean vortex 
rings. The wave was described to mainly affect the Dean vortices but not the boundary 
layers on the wall. For 5042 < Rec 2 < 5270 a secondary Hopf bifurcation was discovered 
leading to a quasi-periodic flow. The second wave, created for Re > Rec 2 , arises as an 
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array of oblique vortices localized in the two boundary layers of the basic flow at the edge 
of the Dean vortex regions. The oscillatory perturbation flows in both the periodic and 
the quasi-periodic regimes were anti- symmetric with re s pect t o the equatorial midplane 
of the torus for <5 = 0.3. For ^ = 0.1 IPiazza fc Ciofalol (|201D found a direct transition 
from a steady to a quasi-periodic flow between 5139 < Re < 5208, associated with hys¬ 
teresis. The traveling waves for (5 = 0.1 were symmetric with respect to the equatorial 
midpla ne of the torus. _ _ 

Like Cioncolini &: Santini ( 2006ll . most studies on curved-pipe flow aimed at reliable 
pressure-drop correlations, because the friction factor is an important quantity for in¬ 
dustrial design. Moreover, all experimental studies on transition to turbulence in fully 
developed curved pipes have been cond ucted in helical pipes w ith a small pitch. While 
the closed flow in a torus considered by Piazza &: Ciofalol ( 2011 ) offers a perfect geome¬ 
try for theoretical and numerical investigations, experiments using a torus are difficult, 
because the driving pressure g r adient cannot as easily be im pose d as in open systems . 
Apart from Madden fc Mullin ( 1994) . del Pino et a,l. ( 20081 ) and Hewitt et a,l. ( 2011 ). 
who experimentally studied the transient flow in a torus during spin-up from rest and 
spin-down from solid-body rotation, we are not aware of any other experimental work on 
the flow in a torus. 

In the present work we experimentally investigate the flow and its instability in a 
torus by means of visual observations, high-speed stereoscopic particle image velocime- 
try (SPIV) and pressure drop measurements. Particular attention is paid to the detection 
and investigation of traveling waves. To that end an experiment was set up where the 
flow inside a torus is driven by a moving sphere. For a sufficiently large aspect ratio 
(circumference-to-pipe-diameter ratio) the perturbations induced by the sphere remain 
localized and bulk flow properties can be measured which are independent of the par¬ 
ticular driving. We measure the full three-dimensional and time-resolved structure of 
the transitional flow including all three velocity-components over the entire length and 
cross-section of the pipe. 


2. Experimental setup 

2.1. The facility 

The toroidal pipe is realized as a toroidal cavity (torus) in a stationary block of plexiglas. 
A ferromagnetic sphere, which is actuated by a strong magnet from the outside, is put 
into the toroidal cavity to drive the flow in the cavity. The main components and the 
driving mechanism of the experiment are sketched in figure [TJ The stationary block of 
plexiglas is built of two disks of highly transparent and polished perspex into which 
a concentric notch of semi-circular cross sections had been machined. The two disks 
are mounted mirror-symmetrically. The two notches in the disks hence form a closed 
toroidal cavity, i.e. a curved tube which is to be filled with the working fluid. To drive the 
fluid motion a ferromagnetic stainless chromium steel sphere with diameter slightly less 
than the toroidal tube is placed in the toroid. The steel sphere is actuated from outside 
the toroidal cavity using a strong permanent magnet mounted on a rotating boom. To 
achieve a constant and precisely adjustable flow rate in the torus the boom is rotated at 
a constant angular velocity, thereby steadily moving the sphere inside the torus, driving 
the fluid motion. The driving system to rotate the shaft consists of an electric gear motor 
combined with a belt drive (both not shown in figure [T]). The Reynolds number is defined 
as Re = Ud/v, where d is the diameter of the tube, v the kinematic viscosity of the fluid. 
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Figure 1. Assembly drawing (side view) of the experimental setup. Lower (1) and upper (2) 
plexiglass disk are bolted (3) and sealed using a 2.5 mm rubber O-ring (4) in a 2.2 x 2.5 mm 
notch. A total of 32 (16 outside and 16 within the torus-notch, in an angular distance of vr/8) 
screws is used to ensure leakproof tightness. A flat pot magnet (5) with threaded stem is incor¬ 
porated to the boom (6) right below the centerline of the tube, leaving an adjustable distance of 
1 mm to the lower plexiglass disk. The boom is rotated around the shaft (7) by a geared direct 
current motor (not shown). The diameter of the torus is D = 614 mm, the diameter of the tube 
is d = 30.3 mm. Drawing not to scale. 


and the bulk velocity U = VlD/2 with angular velocity of the rotating boom and D 
the diameter of the centre circle of the torus. 

When measuring flow details with SPIV a large tube diameter d is desirable. Therefore, 
a tube diameter of d = 30.3 mm was chosen, slightly larger than the diameter of the sphere 
dg = 30 mm, to permit rolling motion of the sphere. This left a small sickle-shaped gap 
of maximum 0.3 mm at the upper half between torus wall and sphere. The total area of 
the gap is Agap = 14.2 mm^ corresponding to 1.9% of the tube’s cross section. Through 
the gap between sphere and tube wall a small leak flow will arise. When the fluid in the 
torus has reached a constant velocity after spin up from rest the sphere has to balance 
just the minor friction losses in the tube which amounted to, e.g. « 20 Pa at Re = 5000. 

Being interested in a small curvature 5 a large diameter D of the torus is preferable. We 
used a centerline diameter oiD = 614 mm, resulting in a curvature S = d/D = 0.049. This 
corresponds to a large aspect ratio (circumference to tube diameter) ofP = t:/ 5—1 ~ (Si. 
To achieve a surface finish of the tube comparable to that of glass pipes the perspex 
surface was polished. The torus dimensions including the manufacturing tolerances are 
D = 614 ± 0.1 mm and d = 30.3 ± 0.03 mm, respectively. 

According to the manufacturer’s specification the pot magnet had ~ 216 N adhesive 
force perpendicular to the surface to which the magnet should stick (equivalent to 22 kg 
bearing capacity). However, the actual adhesive force decreases quickly with increasing 
distance and is also influenced by the alloy of the object of magnetic attraction. Since 
the minimum distance between magnet and sphere in the setup is ~ 5 mm (see figH]), 
the actual adhesive force imposed on the sphere is about 0.49 to 0.98 N (0.05 to 0.1kg) 
according to specification. This was found sufficient to control the sphere in the liquid- 
filled torus. Owing to the high inertia of the rolling sphere compared to any forces caused 
by the flow fluctuations the sphere would always roll very smoothly and accurately follow 
the magnet. Observations of the sphere by a camera mounted on the boom did not allow 
to detect any relative motion between the sphere and the magnet, once the steady rolling 
motion was established. Only in the case of very fast or jerky acceleration or deceleration 
the sphere would not follow the magnet and loose the magnetic bond. 
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To detect the temperature of the fluid in the torus, two small Platinum SMD Flat 
Chip temperature sensors ( Vishay Beyschlag, PTS 0603, 100 fl) are used. Every sensor is 
mounted on the front end of a bracket pin with a diameter of 2 mm. The bracket pins are 
plugged in holes which were drilled into the upper plexiglas disk in an angular distance of 
7r/2, just above the upper apex of the tube. The holes would just not shove through into 
the toroidal volume but leave a thin indentation of approximately 0.08 mm between the 
temperature sensors and the fluid in the torus. For data acquisition the two sensors are 
connected to a Validyne UPC 2100 PCI sensor interface card. The mean value of the two 
sensor signals, averaged over 20 s, is used to determine the temperature of the fluid in the 
torus with an accuracy of ±0.1 K. The ambient temperature in the laboratory could vary 
slowly by approximately ±2 K during a day. As individual measurements would never 
take longer than 5 minutes no additional measures were taken to control the temperature 
of the fluid. To eliminate problems due to the diffe ring refractive indice s such as unwanted 
displacement, hidden regions and multiple images ( Lowe k, KutdHo^ l. arising for objects 
placed inside containers with cylindrical walls, a refractive-index-matc hed mixture of 
distilled water and ammonium thiocyanate was used as working fluid ( BudwigI Il994l : 
Hopkins et al. 2000ll . The kinematic viscosity v of the index-matched fluid as a function 
of temperature was determined by means of a Schott Cannon-Fenske capillary viscometer. 
For temperatures ranging from 20 to 30 °C the kinematic viscosity of the index-matched 
fluid was approximated by the fourth-order polynomial 

1/ = 4.5 X 10"'^ ± 2.7 X 10"^ T - 1.8 X 10“® ± 4.9 x 10"^° - 4.9 x lO"^^ T^, (2.1) 


which was fitted to the measured data by least squares. Here, T is measured in degrees 
Celsius and v in m^/s. A combination of a Lab View-program and a system consisting 
of controller unit (LSC 30/2, 4-Q-DC, maxon motor ag), USB data acquisition card 
(National Instruments 6008) and direct current motor (RE 25, maxon motor ag) with 
speedometer (DCT 22) and planetary gear (GP 32) was used to set and continually con¬ 
trol the Reynolds number. The Lab View-program contained the functional relationship 
between temperature, viscosity and the geometrical data based on which the required 
angular velocity of the sphere for a specified Reynolds number was calculated and imple¬ 
mented. Optionally, an automatic mode can be used to increase or decrease the Reynolds 
number stepwise by predefined step sizes and at predefined times. The setup provided the 
possibility to investigate the flow in the torus for adjustable Reynolds numbers ranging 
from 1 000 to 15 000 with an accuracy of ±2% taking into account the leak flow past the 
sphere ("see 13.1.3]) . 


2.2. Measurements 

For preliminary investigations and visual identification of pertinent flow structures, the 
fluid was seeded with glitter particles (polyester glitter, Sigmund Lindner GmbH) which 
reflect the light and visualize the flow structure. However, the conspicuity was limited, 
because the particles tended to get stuck in the gap between sphere and torus wall and 
would hence block the sphere if their concentration was too high. Therefore, only a limited 
amount of particles could be used. Still photographs as well as video recordings were made 
with a Nikon D 7000 DSLR camera. Movies of the flow structures were recorded with 
a frame size of 1920 x 1080 pixels at 24fps. For the video recordings the camera was 
mounted on an additional boom which rotated at the same speed as the boom driving 
the sphere. This allowed observations in the frame of reference rotating with the angular 
velocity of the sphere which nearly equals the mean flow velocity U. 

The pressure drop between two pressure holes (bore holes with internal diameter of 
1mm) was determined using a differential pressure sensor (Validyne DP103, ultra low 
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Figure 2. Sketch (side view) of the stereoscopic PIV system. To reduce optical reflections, a 
prism is placed on the upper plexiglas disk of the torus. (1) lower plexiglass disk, (2) upper 
plexiglass disk, (3) screws, (4) prism, (5) light sheet, (6) cameras. A few drops of glycerine were 
used to create a thin film between (2) and (4), enabling seamless contact. Drawing not to scale. 


range wet-wet differential pressure transducer). The pressure transducer was calibrated 
using a Betz manometer and it had a full range of 140 Pa with an accuracy better than 
0.35 Pa. The pressure holes were located at the top of the upper side of the tube. The 
sampling rate was 10 ms. 

A SPIV system from LaVision has been used to measure the three components of 
the velocity vectors over a cross section perpendicular to the flow. It consisted of a 
pulsed Quantronix Darwin Duo laser (diode pumped Nd:YLF laser, wavelength 527 nm, 
60 mJ total energy) and two Phantom VIO high-speed cameras with a full resolution of 
2400 X 1900 px. The temporal resolution was set to 100 Hz. Silver-coated hollow glass 
spheres (S-HGS-10, mean diameter 10 /im, p = 1.4g/cm^, Dantec Dynamics) were used 
as seeding particles. __ 

Previous investigators (see e.g. van Doorne fc WesterweefcoOTl) have demonstrated the 
possibilities and potential of SPIV to investigate and capture the appearance and devel¬ 
opment of transitional flow in a pipe. The same principle is applied in the experimental 
setup used for this work. Figure [2] provides a sketch of the setup which was used for 
stereoscopic PIV, showing the arrangement of plexiglass disks, light sheet and 2 cameras 
in a 45° viewing angle. Due to the 45° viewing angle the light would be refracted at the 
outer wall of the plexiglass plate, which would reduce the effective viewing angle of the 
cameras and degrade the image quality. Therefore, a prism made of perspex is attached 
to the upper plexiglass plate such that the optical axis is perpendicular to air-perspex 
interface. 

The calibration target, needed for the calibration of the SPIV measurements, was cus¬ 
tom made by Die Signmaker GmbH (Gottingen). It consisted of 5 mm-spaced lattice 
of black dots with a diameter of 1mm printed on both sides of a 1.5 mm perspex disk 
with a cross-sectional diameter of 30 mm, i.e. slightly less than the tube diameter. The 
disk was kept in position (congruent with the light sheet of the laser) for the calibra¬ 
tion procedure by a tiny locking screw drilled into the torus wall. After the calibration 
images were shot the screw was removed and the calibration target was released. The 
experimental setup was too sensitive to be disassembled after calibration, as dismantling 
and repeated assembly would have caused small displacements and hence spoiled the 
calibration. Therefore, the disk was just left in the torus and went with the flow. After 
a few turns of the actuator the disk would then be located in front of the sphere, not 
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Figure 3. Salient flow structures in a cross section of the toroidal pipe and notation. A cartesian 
{x, y, z) coordinate system is used in the meridional observation plane at constant toroidal angle 
ip. The poloidal (meridional) angle is a. The a;-axis is directed toward the symmetry axis of the 
torus. The symbols hs and Is denote regions of high-speed and low-speed streamwise velocity, 
respectively. A,B,C and D denote the north-pole (B), south-pole (D), outer equatorial (C) and 
inner equatorial (A) points. 


perturbing or changing the flow far from the sphere in addition to the perturbations 
induced by the actuator itself. 

The evaluation of the 3D-vector fields from the PIV images was performed with com¬ 
mercial PIV-software (DaVis 8, LaVision). The interrogation area was 32 x 32px^ with 
an overlap of 50%. The maximum particle displacement between two frames was ap¬ 
proximately llpx in horizontal and 8px in vertical direction. Subsequently, the acquired 
data were analyzed with Mathworks Matlab R2011b which was used as programming 
environment. 

A fibre-flow 2D-LDV system of Dantec Dynamics was employed for laser-Doppler 
velocimetry measurements. From the LDV data we obtained spectra of the streamwise 
velocity at single points within the cross section of the pipe. The average sampling rate 
was 200 Hz. 


3. Results 

Most of the space-resolved measurements are taken in a meridional cross section of 
the torus. Therefore, we use local Cartesian coordinates {x,y,z) as shown in figure [3l 
where {u,v,w) are the respective Cartesian velocity components. The a:-axis is directed 
radially inward towards the center of the torus. Figure [3] also sketches characteristic 
cross-stream flow structures of the basic steady flow, where we use the same notation as 
Webster &: Humohrevl (|l997r) . The steady basic flow is mirror symmetric with respect to 
the equatorial plane ?/ = 0. In a wide core region around the axis y = Q the cross-stream 
flow is directed radially outward in negative a;-direction from the low-speed core region 
(Is) to the high-speed core region (hs) of the streamwise velocity. The cross-stream flow 
returns radially inward in two symmetrically located wall jets called cross-stream wall 
layers. The thickness of these wall jets grows as they turn inward. Before reaching the 
inner equatorial point C the jets turn sharply in region I and merge to form the core 
flow in the low-speed region directed radially outward. 


3.1. Evaluation of the experimental setup 

We are interested in the toroidal pipe flow driven by a prescribed steady volume flux. Due 
to the particular driving by a moving sphere the ideal flow is perturbed in the vicinity 
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of the sphere. To assess to which extent the flow field in the torus can be considered 
fully developed, i.e. stationary and independent of the driving mechanism, the effect of 
the rolling sphere on the flow has to be evaluated. The perturbations arise in form of an 
entrance-length effect and a leakage through the gap between the sphere and the torus. 
As the curvature is constant throughout this work, the strength of the flow is measured 
using the Reynolds number instead of the Dean number. 

3.1.1. Entrance-length effect 

If the sphere moved as a plug with constant velocity without rolling a pure streamwise 
and steady velocity field would be enforced at the surface of the sphere. Due to the 
imposed velocity on the semi-spherical end wall and the velocity discontinuity along the 
line of contact, there exists a return flow in the frame of reference moving with the sphere 
which represents a perturbation to the fully developed toroidal-pipe flow decaying away 
from the sphere. In operation the sphere is rolling in the torus. While the rolling motion 
prevents any perturbation of stick-slip type and thus adds to keeping the translational 
velocity constant, provided the rotation rate of the beam is constant, the rolling motion 
induces an additional cross-stream velocity component which will also decay away from 
the sphere under subcritical conditions. 

To quantify the decay of the above perturbations we measured the streamwise velocity 
w as function of the distance from the sphere using SPIV. Figure S] shows the result 
as function of time (related to the distance from the sphere by U) for Re = 3600, 
Re = 4300 and Re = 4700 in the center of region I at the point (x, y) = (0.17 d, 0.31 d). 
This point has been selected, because the flow instability can be measured sensitively 
in this region. Shown is a period of 12 s, equivalent to approximately one full revolution 
of the sphere around the torus. At Re = 3600 (figure the flow is fully developed, 
except for the immediate vicinity of the sphere where the streamwise velocity component 
varies strongly. In some distance from the sphere the measured velocity w is constant up 
to small fluctuations which are mainly due to the noise of the PIV data. Excluding the 
range disturbed by the sphere, the RMS value of the streamwise velocity is 0.78% of the 
mean value. 

At Re = 4300 (figure |4 }d) a very distinct sinusoidal modulation of the streamwise 
velocity is observed which differs markedly from the signature of the wake flow. With 
the center of the sphere at (/? = 0 the oscillation amplitude is nearly constant in a certain 
region around tp = tt. To analyze the oscillatory flow in the bulk we consider only a time 
interval of 2.25 s corresponding to 0.5 m or ~ 7r/2 which leaves a sufficient safety margin 
from the region disturbed by the sphere. This range is marked as analyzed in figure IHd 
and it is typically located on the side of the torus opposite to the sphere. At Re = 4700 
a dominant frequency can still be recognized, but the amplitude is no longer constant. 

To confirm that the perturbations induced by the rolling sphere are spatially restricted 
we consider the kinetic energy. Decomposing the velocity field into a temporal mean and a 
fluctuation part it = u-\-u' we define the streamwise and cross-stream fluctuation energies 
E|| = {w'^) and E± = {u'^ v'^), respectively, where the brackets indicate averaging over 

the (x, ?/)-plane. Figure [5] shows that the range of perturbed flow is confined to the vicinity 
of the sphere as indicated by the arrows, while further away from the sphere the kinetic 
energy of the fluctuations is nearly constant and very small. For Re = 3600 the flow is 
perturbed for Ri 2.4 s during passage of the sphere, corresponding to an azimuthal angle 
of 0.46 TT slightly asymmetric about the location of the sphere. It is concluded that there 
exists a range of fully developed flow which is nearly unperturbed by the details of the 
flow around the rolling sphere. 

As the Reynolds number increases the unperturbed range decreases. Based on an 
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Figure 4. Time traces of the streamwise velocity w at {x,y) = (0.17 d, 0.31 d) for Re = 3600, 
Re = 4300 and Re = 4700 as function of time t (bottom axis) and azimuthal length 2 (top 
axis). A typical PIV-run of 12 s, comprising approximately one full revolution of the sphere, is 
displayed. The coordinate origin was defined by the transit of the sphere. The vertical dotted 
lines are drawn with a spacing Ap = 7r/2. The velocities of the sphere are Us = 0.19 (a), 0.22 
(b) and 0.24 [m/s] (c). 
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Figure 5. Transit of the sphere through the measurement plane for Re = 3600. The spatially 
averaged (over x and y) kinetic energy of the velocity fluctuations is displayed the direct vicinity 
of the sphere. Shown are the streamwise (F-jlin), the cross-stream (Ui^n) the total fluctuation 
energies Fikin = -f 


extrapolation of the measured length of the zone disturbed by the sphere at low and 
moderate Reynolds numbers the length of the zone of the nearly fully-developed bulk 
flow for Re = 6000 is at least 15% of one full revolution of the torus. This range is 
predicted at G [O.SStt, I.IStt] relative to the sphere at (/? = 0. Even though the lengths 
of the wakes upstream and downstream from the sphere differ, we selected a symmetric 
region for simplicity, leaving a sufficient safety margin to the perturbed regions. All data 
analyzed were taken from these unperturbed angular sections. 

3.1.2. Plunger with sphere 

To assess the effect of the rolling of the sphere on the end-wall-induced perturbations 
and possibly on the bulk flow the driving was modified. Instead of a sole sphere a plunger 
was used made from polyoxymethylene (POM) (Figure [6]). The plunger was machined 
in form of a slightly modihed cylindrical block with a diameter of 30.1 mm by which the 
total area of the gap was reduced to 9.49 mm^ or 1.3% of the tube’s cross section. 

A bore of 20 mm diameter was drilled into the plunger, providing space for a steel 
sphere with 18 mm diameter. This steel sphere was also controlled by the magnet from 
outside, but its rolling motion would not influence the flow in the torus. As the plunger 
was cylindrical, its shape did not perfectly match the slightly curved torus. Due to the 
lubrication of the fluid between the two solid surfaces of tube and plunger, the plunger 
would slide through the tube very well without excessive stresses or seizures at the wall. 
Only very mild wear, mostly originating from the tracer particles in the fluid, was noticed. 

We found that the design type of the actuator had very little influence on the length 
of the perturbation flow near the moving end wall and no influence on the bulk flow. For 
that reason only a rolling sphere was used for all measurements. 


3.1.3. Leak flow past the sphere 

The azimuthal velocity Us of the center of the sphere is precisely known as it is con¬ 
trolled by the angular velocity D of the rotating boom. Comparing the associated flow 
rate Ugr^df /A with the measured flow rate allows to determine the leak flow past the 
sphere. The actual mean velocity in the torus U was determined in two ways. First, the 
pressure drop Ap over the angle Ap = 7r/2 was measured . The result was compared t o 
pressure-drop data obta ined from published correlations of White! 1 1929ll : Hasson ( 1955 1: 
Mishra fc Guptal ( 1979ll for the laminar regime, which all deviate less than 1.5%. Up to 
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Figure 6. Sketch of the first type of plunger within the torus. A small sphere of ds = 18mm 
placed inside the plunger is used as object of magnetic attraction to move it. Only the streamwise 
moving walls of the plunger, instead of the surface of the rolling sphere, interact with the flow 
in the torus. 

Re < 2500 we find excellent qualitative agreement with these correlations but a practi¬ 
cally constant deviation which can be attributed to a slightly slower velocity of the fluid 
U than the azimuthal velocity Us of the sphere by ~ 3%. In addition, U was obtained 
from SPIV-measurements. The resulting leak flux amounted to 3% of the actual mean 
volume flux, consistent with the flow rate from the pressure drop. Therefore, we can 
safely determine the mean velocity as U = 0.97 C/s. It must be noted, however, that a 
slightly non-linear behavior of the leak flow for increasing Reynolds number could not 
be ruled out. 


3.2. Velocity-field measurements 

3.2.1. Steady basic flow 

In the following the flow field in the bulk is considered which is practically unaffected 
by the end-wall perturbations. For small Reynolds numbers the flow field is unique and 
steady. For Re = 3600 the steady basic flow has developed its characteristic structure 
which is caused by the pipe curvature. The velocity field measured by SPIV in the 
entire cross-section is shown in figure 0 As the fluid in the core region is centrifugally 
driven toward the outer wall, the maximum streamwise velocity is located in the plane 
of symmetry and near the outer wall. Therefore, the gradient of the streamwise velocity 
component U|| is very high near the outer wall. As a result of the secondary, radial outward 
motion in a wide range about the symmetry plane and due to continuity, fluid with high 
streamwise momentum is transported radially inward along very thin cross-stream wall 
layers. This is visible from the cross-stream velocity component u± and the elongated 
contour shapes of the streamwise velocity. As the wall jets evolve radially inward they 
widen. At approximately a = ±7r/4 the jets turn and merge with the broad radially 
outward stream, without separating from the wall. The two symmetrically located turning 
regions of the secondary flow, denoted region I, represent the vortex cores of the two Dean 
vortices. 

The stationary waviness of the radial outward cross-stream flow near the equatorial line 
in figure [7] indicates a weak symmetry breaking. It could not be determined conclusively 
whether the small loss of symmetry is due to the flow physics or due to systematic 
measurement errors. Possible sources of error are the symmetry breaking by the rolling 
of the sphere and an imperfect calibration due to the custom made calibration target (see 
sec. IQ) . However, since the cross-stream velocity is about one order of magnitude smaller 
than the streamwise velocity it is evident that the deviation from perfect symmetry is 
diminutive. 

Figure [8] shows profiles along the vertical a; = 0 of the streamwise velocity u,|| = w{y) 
and the absolute value of the cross-stream velocity u±_ = \/u^{y) + u^(y) (~ \u{y)\ for 
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x/d 


Figure 7. Entire cross-sectional flow fleld of the tube in the steady laminar regime (at 
Re = 3600). Streamwise (color contours) and in-plane velocity (vectors), both normalized with 
the mean streamwise velocity U. To eliminate measurement noise, the flow field is an average 
of 400 PIV-images. 


a; = 0) for Re = 2400. The wall layers are clearly visible with global maxima of the 
absolute value of the cross-stream velocity at y = ±0.46 d in a distance of 0.04 d from 
each wall. The two other maxima of the absolute value of the cross-stream velocity at 
y = ±0.34 d mark the edge of the cross-stream return flow in the interior. Between these 
extrema of u± the streamwise velocity rt|| reaches its extrema at y = ±0.38 d before 
dropping sharply towards the wall. The boundary layers become thinner with increas ing 
Reynolds number and/or curvature of the pipe (see also lWebster fc HumphrevlIlQQTfl . 


3.2.2. Mean-velocity profiles for increasing Reynolds number 

For higher Reynolds numbers the basic flow becomes unstable and time-dependent. 
To characterize the flow in a wide range of Reynolds numbers we consider the time- 
averaged streamwise velocity ?1|| = w/U in the two orthogonal planes x = 0 and y = 0. 
Figure ini shows the profiles U||(a;,0) (denoted by h for horizontal) and i2||(0,2/) (denoted 
by V for vertical) for various Reynolds numbers. Characteristic for curved-pipe flow are 
the flat vertical profiles and the nearly linear variation of the horizontal profiles. As 
the Reynolds number increases the boundary layers become thinner, associated with 
increasing strain rates at the walls. Moreover, the slope of the horizontal velocity profiles 
becomes smaller, associated with a reduction of the velocity extrema which is most 
significant for Re = O(IO^). It is remarkable that the horizontal profiles of the mean 
velocity do not change very much between Re = 3600 and Re = 6000, even though the 
flow for Re = 6000 is strongly time-dependent. In addition to the profiles of the mean 
streamwise velocity along y = 0 and a; = 0, U||(0.25 d, y) (denoted by v*) is shown. As 
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Figure 8. Vertical {x = 0) profiles of the streamwise (u||) and in-plane (u±) velocity 
magnitudes for Jie = 2400. Velocities not to scale. 

can be seen, the mean velocity profiles along x = 0.25 d (v*) do experience a change 
between Re = 3600 and Re = 6000. This change is related to the the onset of oscillatory 
flow whose amplitude is largest near the maxima of along x = 0.25 d. 

3.3. Primary instability: The onset of time-dependence 

When the Reynolds number is increased the steady basic flow becomes unstable and 
time-dependence sets in at Rec = 4075. The supercritical flow arises as a wave which 
travels downstream with a phase velocity which is slightly larger than the mean flow. 
With increasing Reynolds number the amplitude grows continuously from zero at the 
threshold. Figure fTOl shows a short section {Aip = 0.137r) of the torus at Re = 4350. The 
wavelength A is indicated by arrows. The wave can be recognized much better in the 
short movie mov-RefSBO.ltum.avi (supporting material, online) than in the image. 

From visual observations the wavelength A of the finite-amplitude wave in terms of 
the arclength along the center of the pipe can be estimated as A = A(/3iA/2 k, (0.06- 
0.08)7r£)/2 Ki (2-2.5) d for Re A 4350. In the range of 4075 ^ Re ^ 4350 the wave 
is almost stationary in a frame of reference moving with the sphere. The wave celerity 
c = (1.1-1.13) U is about 10% higher than the mean velocity. For Re ^ 4350 the strict 
periodicity of the wave is lost and the peak-to-peak distance of the disturbed wave is 
varying between 2-4 d. Moreover, for Reynolds numbers Re ^ 4580 the wavy motion is 
interrupted by short irregular bursts. For even higher Reynolds numbers. Re ^ 5050, no 
dominant period can be detected by visual observation. 

At Re = 4300 the finite-amplitude wave is already well established. Figure [TT] shows 
the streamwise velocity w in the (x, 2 /)-plane at four instants during one period. Each 
image is obtained from averaging eight successive periods. The oscillation amplitude of 
the cross-stream velocity is very small (not shown). The fully time-resolved oscillatory 
flow is shown in Re4300.timeresvelfield.gif (supporting material, online). The oscillations 
are most pronounced in the vicinity of the contour level 1 U, in particular in region I (see 
figure [SD. During the oscillation the finger of high streamwise momentum in region I is 
oscillating back and forth in meridional direction (parallel to the wall). When the fingers 
elongate towards point A at the inner wall, the radially stratified streamwise momentum 
in the center of the pipe moves towards point C on the outer wall, the shift however 
being smaller, and vice versa. 
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Figure 9. Horizontal (li) and vertical (v) profiles of the mean streamwise velocity at j/ = 0 and 
X = 0, respectively, for different Reynolds numbers, v* denotes vertical profiles along x = 0.25 d. 



Figure 10. Still picture from a movie (see mov-Re4350.1turn.avi) of a short section of the torus 
from above at Re = 4350. The camera was following the flow with the velocity of the sphere Ua ■ 

Figure [12] shows the vector field of the cross-stream velocity fluctuations and selected 
contour lines of the streamwise velocity fluctuations at eight instants during one period 
(see also Re4300.fluctvelocity.gif in the supporting material, online). Each image is ob¬ 
tained from averaging 8 periods of oscillation. Despite of the slight asymmetry of the basic 
flow the oscillations (deviations from the temporal mean) are symmetric with respect to 
the horizontal plane. The magnitudes of the velocity oscillations in the streamwise and 
the cross-stream directions are of the order of ±0.061/ and ±0.021/, respectively. The 
cross-stream velocity field in (a) consists of a radial outward flow from the center (la¬ 
beled 5) of the pipe in the right half of the cross section. The radially outward flow at 
a ft! 7r/4 towards and across region I is particularly strong (point 2). Near point 6 at 
the edge of the cross-stream wall jet at a ft 37r/4 the cross-stream velocity is oriented 
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Figure 11. Contours of the (instantaneous) streamwise velocity w over one wavelength A at 
Re = 4300 in increments of 0.05 [/. a) A = 0, b) A = 1/4, c) A = 1/2 and d) A = 3/4. The level 
of 1.55 U and 117 is indicated for reference. 

tangentially in negative a-direction. The cross-stream oscillation amplitude is very small 
in the regions of low and high streamwise basic-state velocity near the equatorial points 
A and C, respectively. Half a period later (e) the cross-stream-flow field is reversed. These 
cross-stream velocity perturbations act on the underlying mean streamwise velocity field 
(basic flow). Owing to the cross-stream gradients of the mean streamwise flow the cross¬ 
stream perturbation flow creates the streamwise perturbations (labeled by numbers in 
figure [T^ which can be considered as streaks. The periodic cross-stream velocity acts, 
in particular, on the fingers of the streamwise mean flow. Since the streamwise mean 
velocity exhibits a minimum and a maximum along the ray a = 7r/4 (see figure [TT]) the 
radially outward cross-stream perturbation flow along this direction creates three streaks 
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Figure 12. Fluctuating velocities over one wavelength \ aX Re = 4300. Vector field of the 
cross-stream velocity fluctuations and contours of the streamwise velocity fluctuations, both 
normalized with the bulk velocity U. For the sake of clearer representation w' is displayed only 
for ±0.008 to ±0.03 in increments of 0.004. The maximum levels would be ±0.065. 
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(a) View from the inside of the bend 



(b) View from the outside of the bend 



Figure 13. Isosurfaces of streamwise velocity oscillations (streaks) at ±0.0281/ for Re — 4300 
shown over one wavelength A. Two views from opposing sides (I denotes inner wall, O outer 
wall) are shown. The flow is from left to right. Blue indicates a negative, light gray indicates a 
positive streamwise velocity perturbation. The images have been averaged of eight periods. For 
reasons of representation, the velocity field is mapped to a straight pipe. 


(1,2,3) of alternating sign. During the temporal (or spatial) evolution, the central streak 
(2) grows while the other streaks move radially inward (3) and outward (1). After streak 
(2) has grown to a considerable size (figure[T2ji) it splits and the original streak structure 
(figure fT^ ) is recovered in figure IT^ . albeit with a different sign. The pair of radially 
inward moving streaks merge at the equatorial plane (figure [T2b"c) and subsequently 
merge with the two merged streaks (6) which have also moved towards the equatorial 
plane. 

The time- and space-resolved structure of the streamwise velocity fluctuation is shown 
in figure [13] for Re = 4300. Isosurfaces at ±0.02817 are shown from the inside (a) and 
from the outside of the bend (b). The streaks, i.e. the regions where the streamwise 
velocity is lower or higher than the steady laminar value, are shown in blue and light 
gray, respectively. The mirror symmetry with respect to the equatorial plane is obvious. 

The spatial distribution of the fundamental Fourier mode Ai(x,y) with frequency /i 
is shown in figure [T4l for the streamwise velocity fluctuations U||. The fundamental mode 
is dominating at Re = 4300 and the corresponding streamwise velocity fluctuation is 
sharply peaked at (x, y) = (0.17 d, ±0.31 d). Near the center another region of streamwise 
velocity perturbation is found which is quite wide, but with smaller amplitude. The 
streamwise velocity fluctuation is practically absent in the low-speed core region and 
in the cross-stream wall layers. This result underlines the importance of r egion I for the 
instability, consistent with the observation of lWebster fc Humphrey ( 19971) who, likewise, 
found a maximum oscillation amplitude in this region of the flow. 

The absolute maximum of the velocity perturbation in region I provides the best signal- 
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Figure 14. Spatial distribution of the amplitude of the fundamental Fourier mode A\{x,y) 
in the streamwise velocity-fluctuation field for Re = 4300 measured by PIV. The amplitude is 
given in arbitrary units. 


to-noise ratio for measuring the streamwise velocity. Since the location of the maximum 
was nearly independent of the Reynolds number, the position (x^y) = (0.17d, 0.31 c?) 
was selected for LDV measurements. The measured streamwise velocity perturbation 
w was Fourier analyzed for different Reynolds numbers. Figure [15] shows spectra at 
four different Reynolds numbers obtained from the range ip = [O.SStt, I.IStt] ((/j = 0 
corresponds to the location of the sphere) during a single revolution of the sphere. The 
dimensionless frequency (Strouhal number) is obtained as / = /d/C/, where / is the 
frequency in Hz. At Re = 4350 (a) a fundamental frequency fi and its weak second 
harmonic frequency /2 = 2/i is clearly resolved. At Re = 4400 (b) and Re = 4700 (c) 
the same fundamental frequency fi is still dominant, but its amplitude is less than for 
Re = 4350. Furthermore, additional frequencies in the vicinity of /i can be detected. At 
Re = 5000 (d) no additional frequencies apart from fi and its higher harmonics can be 
resolved. 

The spectrum for Re = 4350 (figure [TSb l is reproducible for each revolution of the 
sphere and it is characteristic for the range 4075 ^ Re ^ 4350. The amplitudes wi ^2 differ 
slightly from revolution to revolution. The dominant frequency /i, however, remained 
constant. For cases (b), (c) and (d) the measured spectra varied distinctly from revolution 
to revolution (not shown). 

The first occurrence of fi within the velocity spectrum indicates the first critical 
Reynolds number. The dependence of the amplitude wi on the Reynolds number is 
indicative of supercritical bifurcation. To reduce the deviations among the results ob¬ 
tained from different revolutions, the angular range [0.857r, 1.157r] was Fourier analyzed 
and the result averaged over five revolutions of the sphere for each Reynolds number. 
The result is shown in figure 1161 The error bar indicates the maximum deviation of an 
individual amplitude from the mean value. Up to Re = 4350 the measured data can be 
fitted to a square-root law (black dashed line) 

wi {Re) = ai^J Re — Re^i (3-1) 

with ai = 4.7 x 10“^ and Reel = 4075. The amplitude of the second harmonic W 2 is 
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(a) Re = 4350 (b) Re = 4400 
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(c) Re = 4700 
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Figure 15. Spectra of w measured by LDV at (a;,y) = (0.17d, ±0.31 d) obtained from 
ip = [O.SStt, 1.157r] during a single revolution of the sphere. Shown are the Fourier amplitudes of 
the streamwise fluctuations normalized with the mean velocity in the torus, i.e. w{f)/U. 


also shown. This supercritical bif urcation for 6 = 0 . 049 is in contrast to the subcritical 
bifurcation for <5 = 0.1 found by IPiazza fc Ciofalo ( 201lli numerically in the periodic 
torus. 

Figure ITTI displays /i and /2 in the range 4000 ^ i?e ^ 5400. The fundamental fre¬ 
quency (black squares) decreases slightly with the Reynolds number and can be fitted to 
a linear curve (black dashed line in figure ITT)) 


A = ^ = 1.05 - 1.42 X IQ-'^Re. (3.2) 

Up to the measurement error this behavior is also satisfied by the second harmonic 
frequency / 2 . Apart from /i and /g no further dominant frequencies are found in the 
velocity spectra in the range Reel ^ Re ^ 4350. However, the fundamental frequency 
fi{Re) can be traced to much higher Reynolds numbers into the regime of more complex 
flows in which periodicity in the bulk is lost. 


3.4. Beyond the first instability 

The Hopf instability at Red = 4075 is breaking the continuous rotational invariance 
with respect to p while preserving the mirror symmetry with respect to the equatorial 
plane. The deviation of wi(Re) from the square root law in figure ITHl suggests that a 
secondary bifurcation occurs near Rec 2 ~ 4400. Around this Reynolds number the peak- 
to-peak distance of the flow oscillations was visually observed to start varying between 
approximately 2d and dd, indicating a loss of simple periodicity. The widening of the 
spectrum near fi in figure [THlb and c indicates that a further bifurcation to a state with 
two possibly incommensurate frequencies A and f[ takes place. The second frequency /{ 
is conjectured to exhibit a value below f\. Despite of careful investigation of the velocity 
spectra /( could not be pinpointed as every revolution of the sphere yielded different 
spectra. 
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Figure 16. Normalized amplitude wi/U of the dominant frequency fi (squares) and its har¬ 
monic W 2 IU (dots) for increasing Reynolds number. Error bars represent the maximum devi¬ 
ation of the measured amplitude for single revolutions. The dashed and dotted lines represent 
square-root and linear approximations, respectively. 
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Figure 17. Strouhal number /i (black squares) and its harmonic /2 (red dots) as functions of 
the Reynolds number. Linear fits are indicated as dashed lines. 


Characteristic of the conjectured second instability is a breakdown of the mirror sym¬ 
metry with respect to the equatorial plane. This is conhrmed in figure [18] which displays 
instantaneous flow fields at Re = 4600. The symmetric oscillations seem to ’break down’ 
from time to time and, associated with the break down, the symmetry with respect to 
the equatorial plane is lost temporarily. The asymmetry is observed in the interior as well 
as in the wall layer. The asymmetry in the interior arises in form of a non-zero velocity 
component u(x,0) along the x axis in the outer half of the pipe {x < 0) and visible 
interior vortex structures for x < 0 which are asymmetric with respect to y = 0. 

The equatorial symmetry is lost due to small-scale flow structures in the wall layers 
near the equatorial plane and also by interior-flow structures. The small-scale vortices 
arise spontaneously in the wall layer. One example is shown in figure[T8|a. A low-velocity 
streak is deforming the isolines of the streamwise flow in the wall layer just below the 
x-axis near point C. Such a streak is usually caused by a pair of cross-stream vortices 
(only one of which can be identified in the total flow in figure fTSla). The structure is 
visible for a short period of time of about 0.2 s. Immediately after vanishing the same 
structure appears for the same short period of time in mirror symmetrical location in the 
upper half of the torus ffigure flSlb). In addition to these localized structures which seem 
to arise erratically, we find asymmetric flows in the interior which seem to be related 
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(a) (b) (c) 
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Figure 18. Cross sectional flow field for Re = 4600 at different moments of time (same as in 
fig. I19II : t = 0.60 (a), t = 0.88 (b) and t = 3.23 (c). Note the small-scale vortex in the wall layer 
as well as the flow asymmetry in the interior. 
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Figure 19. Time traces for Re = 4600 of the streamwise velocity during four seconds (corre¬ 
sponding to « 50% of one full revolution of the sphere or the range p € [7r/2, 37r/2]. Signals are 
shown for five locations PI to P5 as explained in the text. 

to the underlying symmetric oscillations with frequency fi. From time to time interior 
vortices are created in the cross-stream flow in the outer half of the torus {x < 0). An 
example is shown in figure [T 8 ]c. The asymmetric interior-flow structures do not seem 
to be related to the small-scale vortices in the wall layer, but both can be observed for 
Re > 4400. Owing to the limited observation time we were not able to exactly pinpoint 
a critical Reynolds number Rec 2 for the symmetry breaking. 

In order to further characterize the temporal and spatial behavior of the flow beyond 
the second instability figure [19] displays the streamwise velocity w{x, y, f) during a period 
of four seconds for Re = 4600 corresponding to the azimuthal range ip G [ 7 r/ 2 , 37 r/ 2 ]. Sig¬ 
nals shown are obtained from three locations in the mid-plane: PI: (a;i, yi) = (—0.38 d, 0), 
P2: {x 2 ,y 2 ) = (—0.25d, 0), P3: (xs^ys) = (—0.09d, 0), and from two symmetrical loca¬ 
tions in region I: P4: ( 2 : 4 , 2 / 4 ) = (0.17d,0.31 d) and P5: ( 2 : 5 , 1 / 5 ) = (0.17d,—0.31 d). For 
the periodic oscillatory flow for Red < Re < Rec 2 all these signals are periodic as in 
figure m Moreover, the signal from the antisymmetric points P4 and P5 would be iden¬ 
tical. For Re > Rec 2 as shown in figure [T9| we observe the loss of symmetry indicated 
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(a) Re = 6000 
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(b) Re = 9000 
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Figure 20. Time traces of the streamwise velocity during 3 seconds (corresponding to ~ 50% 
of one full revolution for (a) and ~ 70% for (b) respectively) at Re — 6000 (a) and Re = 9000 
(b) for 5 different points within the cross section. PI — P5 as in figure [TO] 


by the differences between signals from P4 and P5. While the flow is symmetric and 
periodic for a considerable time, it is interrupted by seemingly random anti-symmetric 
and high-intensity bursts. Such bursts can be clearly identified in the left half of the torus 
by the signals of PI, P2 and P3 at t = 1.5 and of P2 and P3 at t = 3.3, leading to an 
asymmetric flow as indicated by the differences in the signals from P4 and P5. 

Figure [20] shows the five signals at the same monitoring points for Re = 6000 (a) and 
for Re = 9000 (b). For Re = 6000 the time traces of the streamwise velocity are similar to 
those at Re = 4600. The symmetry of the signals P4 and P5 is only slightly perturbed. 
However, the nonlinearity in the oscillations is stronger which also shows up by higher 
Fourier components in the spectrum (not shown). Furthermore, periodic oscillations are 
difficult to be recognized in the signals from PI, P2 and P3. 

At Re = 9000 (figurel^Db) the flow has changed significantly as compared to Re = 6000. 
The signal is dominated by high-frequency and random oscillations in all monitoring 
points. The symmetry between signals P4 and P5 is completely lost. The character of 
the fluctuations evidence turbulent characteristics. It is noted that the overall amplitudes 
of the fluctuations have decreased at P4 and P5 as compared to lower Reynolds numbers. 

A comparison between figures [20h and b shows that the high-frequency fluctuation 
in (b) are essentially absent in (a). Therefore, we consider the flow at Re = 6000 to be 
chaotic while the flow at Re = 9000 is considered turbulent. For a further classification 
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Figure 21. Friction factor measured in the closed torus as a function of the Reynolds number. 
Straight lines indicate laminar- and turbulent-flow correlati ons fcL a nd fcT, respectively, of 
iMishra &: Gupts] 11197911 and the turbulent-flow correlation of lltd ll 19591 1. both for helical pipes. 


of different flow regimes a more detailed investigation of the turbulent fluctuations is 
required. 


3.5. Friction factor 

The mean streamwise pressure gradient in the range 1000 ^ Re ^ 15.000 was obtained 
by measuring pressure differences between two pressure holes in an angular distance 
Aip = 7r/4. Measurements were only made while the sphere was moving on the side 
of the torus opposite of the bore holes. For Re < 18.000 the pressure difference Ap = 
p(7r/2 -I- tt/8) — p(7r/2 — tt/S) was found to be constant during a sufficiently long time 
span, i.e. any pressure fluctuations were sufficiently small to measure the mean pressure 
difference Ap with a maximum r.m.s. value of 0.4 Pa. 

Figure [21] shows the Fanning friction factor 


^ 2pC/2A(p’ 


(3.3) 


determined from pressure-loss measurements (squares). The measurements are in good 
agreement with published data for helical pipes, in particular, for the laminar regime. 
Deviations from the la minar friction factor fcL approximated by the power law of 
Mishra fc: Gupta ( 1979ll become appreciable for Re > 2000 beyond which the measured 


friction factor decays more rapidly than the power law. This trend is alleviated by the 
onset of oscillations at Red at which point the slope of the friction-factor should change. 
However, we have not enough data to unambiguously resolve the change of the slope. 
After a transitional range a further change of the slope of / can be clearly recognized at 
Re ~ 8000. For Re > 8000 th e measured data are in reasonable good agreement with the 
correlation fcT of Ito ( 19591) for helical pipes. 


4. Discussion 

The torus experiment revealed a first instability of the toroidal-pipe flow in the bulk 
to a traveling wave at Re = 4075. As the Reynolds number is increased from its critical 
value the amplitude of the fundamental Fourier mode grows continuously and homoge¬ 
neously, indicating a forward Hopf bifurcation. For Re > 4400 the fundamental frequency 
develops sidebands. We anticipate that the sidebands are caused by an incommensurate 
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frequency. The traveling wave represents a large-scale structure occupying the full cross 
section of the pipe. In addition, we find intermittent small-scale and nearly streamwise 
vortices which are localized in the wall layer in the vicinity of the impingement point C. 
For Re > 8000 the spectrum becomes broadband and the signal exhibits high-frequency 
fluctuations. The corresponding threshold Reynolds number bracketed by PIV was con¬ 
firmed by pressure-drop measurements which revealed a continuous dependence of the 
pressure drop on Re with a sudden reduction of the slope at Re ~ 8000. 


Sinc e flow instabilities in toroidal pipes have merely been considered bv iPiazza fc Ciofalo 


(|201lh . albeit for periodic bounda ry conditions, our results can onl y be further compared 


with the resu lts for helical pipes of Sreenivasan &: Strvkowski ( 1983ll and IWebster fc: Humphrey 


( 199,81 I 997 I I. Regarding the frictio n factor measuremen t s in th e transitional regime we 


compare our results with those of ICioncolini fc Santinil (j2006f ) for helical pipes. When 


comparing with helically coiled pipes it is tacitly assumed that the particular driving 
mechanism in the present experiment does not sig nificantly modify the bulk flow, nor 
does the small pitch of the helical winding (see e.g. iBerger et al.l[l983t Piazza fc Ciofalol 


201111 . Common to all these investigations, the loss of stability of the steady basic flow 


occurs in the range 4000 ^ Re^i < 6000 with a subsequent traveling-wave regime for 



and quasiperiodic low frequency oscillations near the inner wall of the pipe and high- 
frequency bursts near its outer wall. The signals were detected by two hot wire probes in 
the midplane and at 0.125d from the wall at points A and C. The slightly la rger critical 
Reynolds number of Re = 4200 obtained bv Sreenivasan fc Strvkowsk i 1983ll may be due 


to the small difference in 5 and/or the small amplitude of the traveling-wave at the loci 
of the hot-wire probes (cf. hgure 11411 . Furthermore, the low-frequency signals measured 
bv ISreenivasan fc Strvkowskil l 198,811 at Re = 5000 have neither a constant amplitude nor 


a constant frequency. In fact, they are very similar to our result for the quasiperiodic 
regime (see e.g. P5 in figure uni- Also the coexistence of a low-frequency modulated 
wave with short, high-frequency bursts in different regions of the flow for Re = 5870 
is consistent with the present findings. As Sreenivasan fc Strvkowski ( 1983ll did neither 


s pecify any wavelengths nor am plitudes a more quantitative comparison cannot be made. 

IWebster fc Humphrey (1993) investigated a nominally fully developed flow through a 
helically coiled pipe with a curvature of 5 = 18.2“^ = 0.055 using LDV. The curvature 
is even closer to the present value. In the range 5060 ^ Re ^ 6330 the authors found 
periodic flow oscillations with a constant fundamental Strouhal number / = 0.25 mea¬ 
sured in the inner half of the pipe cross section. This is qualitatively compatible with 
our result. However, we find a nearly linear variation of the fundamental frequency from 
f{Re = 5060) = 0.32 to f{Re = 6330) = 0.15. 


Wh en comparing with the above investigations it must be noted that lSreenivasan fc Strvkowski 
used a long straight -pipe section before the inlet of the coiled pipe, whereas 
Webster fc Humohrevl (jl9931 1 used a ’flow straightener’ in the straight-pipe section, about 
thirty pipe diameters upstream of the inlet of the curved pipe. For Reynolds numbers at 
which the flow is turbulent in the straight pipe two experiments may represent quite dif¬ 
ferent inlet conditions for the curved pipe. Th is could be the reason why th e first critical 
Reynolds number Red = 5060 determined by Webster fc Humphrey! (Il993ll is well aboye 
the critical Reynolds number determined by ISreenivasan fc Strykowskil ( 1983 1. and also 
larger than the present critical Reynolds number Red = 40 75. 

The numerical simulations of lWebster fc Humphrey! ( 1997 1 revealed that the amplitude 






































































Experimental investigation of transitional flow in a toroidal pipe 


25 


of the traveling wave is very small near the midplane and that the maximum stream- 
wise velocity perturbations arise in the center of region I . This is confirmed by our 
measurements (compare figure 7 of Webster fc Humphrey ( 1997f l with figures [T^ and 
m. Wavelengths and phase velocities are of the same order of magnitude, even though 
differences remain. 

The only study to date of the transition to turbulence in a torus is due to Piazza fc Ciofalol 
( 201lh . For curvatures S = 0.1 and S = 0.3 they numerically simulated the flow in a torus, 
driven by an artificial azimuthal body force, for periodic boundary conditions and several 
Reynold s numbers in the range from 3500 to 14700. Upon an increase of the Reynolds 
number IPiazza Sz Ciofalol (20111 found stationary, periodic, quasi-periodic and chaotic 
flows. The flow states found are in qualitative agreement with our experimental results. 
However, there are differences. 

For S = 0.3, a value much larger than the present curvature S = 0.049, Piazza fc Ciofalol 
( 201lh reported a supercritical Hopf bifurcation at Rec = 4575, giving rise to a traveling 
wave along the Dean vortices. A second Hopf bifurcation between 5042 < Re < 5270 
led to a quasi-periodic flow where the additional frequency is associated with oblique 
vortices located near the outer equator at C. While the sequence of bifurcations and 
the regions where the two waves have a sizable amplitude agree qualitatively with our 
res ults, the symmetry with respect to the midplane differs, because both modes found 
by iPiazza fc Ciofalol (|201lll are anti-symmetric with respect to the equatorial plane. 

For 5 = 0.1, a value much closer but yet a factor of two larger than the pre sent 
curvature, the periodic and quasi-periodic waves found bv IPiazza fc Ciofalol (20111 are 
mirror symmetric as in our measurements, but the first instability at Rec = 5175 is a 
saddle-node bifurcation associated with a subcritical Hopf bifurcation. In contrast, we 
find a supercritical Hopf bifurcation for 6 = 0.049. This difference might be attributed to 
the different curvatures investigated. It is also not precluded that the different boundary 
conditions lead to different finite-amplitude flow patterns when different modes of the 
corresponding linear-stability problem are only weakly damped. 

IPiazza fc Ciofalol (201 lll have sketched a tentative map of flow states in the (<5, Re)- 
plane (their figure 30). Accordingly, the range of chaotic flow is strongly stabilized as S is 
increased from zero, giving way to symmetry-breaking bifurcations at larger S. In order 
to elaborate their flow map, further experimental investigations for smaller S would be 
desirable. 


Despite of the differences in the setup and the boundary conditions of the above studies 
on curved pipe flow, there exists no sharp jump in the friction factor upon transition to 
turbulence as is known from the flow in straight pipes. In the present toroidal pipe flow 
the flow changes gradually and the friction factor is a monotonous function of Re. Thus 
no distinct onset of turbulence can be detected by pressure drop measurements. However, 
we found indications for an onset of chaotic motion at Rec 2 ~ 4400 and a clear change 
of the slope of / in the logarithmic friction-factor diagram at Re ~ 8000 which we could 
correlate, by PIV measurements, to the onset of irregular high-frequency fluctuations. 


For helically coiled pipes Cioncolini fc Santinil ( 2006l l analyzed the influence of coil 
curvature by examination of friction factor profiles. Their findings are in very good qual¬ 
itative agreement, as they report a gradual discontinuity in the friction factor profiles, 
marking the end of the turbulence emergence process. It must be noted, however, that 
their results indicate the discontinuity at Re « 7300 for 6 = 0.049 (according to their 
equation (13)), which is considerably lower than Re ps 8000 found in the present inves¬ 
tigation. 
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5. Summary 


A novel experiment has been set up which realizes a flow with a precisely adjustable flow 
rate in a toroidal pipe, allowing to accurately measure the cross-sectional velocity field 
using optical techniques. The facility has proven a useful tool to explore the transition 
to turbulence in a curved pipe. 

Using SPIV distinct flow states have been identified and characterized for the curvature 
6 = 0.049. Mean velocity profiles along the vertical and horizontal axes of the pipe cross 
section have been measured as function of the Reynolds number. Moreover, instantaneous 
flow fields were obtained for different Reynolds numbers. The first critical Reynolds 
number Red = 4075 ± 2% has been determined very accurately. The bifurcation is of 
supercritical Hopf type. A further critical Reynolds number was found at Rec 2 ~ 4400 
at which the flow becomes presumably quasi-periodic. The space resolved measurements 
allowed to identify the spatial and temporal structures of the dominant oscillatory mode 
just above the threshold Red as well as large and small scale fluctuations which arise in 
different regions of the flow for Re > Rec 2 - The SPIV measurements have shown that the 
flow becomes turbulent at Re ~ 8000, a value consistent with the shape of the friction 
factor curve which has been measured in the range of 1000 ^ Re ^ 15.000. 

The present investigation, using a single curvature, has confirmed the scenario of bi¬ 
furcations that precede the onset of turbulent flow. Perhaps most importantly, we have 
measured the traveling wave resulting from the first supercritical instability fully resolved 
in the cross-sectional plane and in time. These measurements allowed to identify the in¬ 
terplay between streamwise vortices and streaks on which the traveling wave is based. 
From these observations we deduce that the first instability is caused by the alternating 
radial gradients (under a ~ 7r/4) of the streamwise velocity near region I. This result is 
in some contrast to previous interpretations in terms of a centrifugal instability. 

In view of the sparsity of detailed results on transitional toroidal pipe flow and the 
qualitative and quantitative differences among the few previous investigations a more 
systematic variation of the curvature would be desirable in order to clarify the remaining 
open points. Of particular interest would be the dependence of the first critical Reynolds 
number on the curvature, the existence ranges of other flow regimes which arise through 
higher-order bifurcations, and the competition with turbulent spots which may arise on 
a decrease of the radius ratio. Another interesting aspect wou ld be the continuation of 
the unstable wave solu tions of the Navier-Stokes equations ( Faisst fc Eckhard^ 2003 : 
Wedin fc Kersw^ 2004 ) as the curvature is increased from zero to curved pipe flow for 


which stable waves are known to exist. 
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